Theory of Inelastic Electron Tunneling from a Localized Spin in 

the Impulsive Approximation 

M. Persson 
The Surface Science Research Centre, 
The University of Liverpool, Liverpool, L69 3BX, UK 

Abstract 

A simple expression for the conductance steps in the inelastic electron tunneling from spin 
excitations in a single magnetic atom adsorbed on a non-magnetic metal surfaces is derived. The 
inelastic coupling between the tunneling electron and the spin is via the exchange coupling and is 
treated in an impulsive approximation using the Tersoff-Hamann approximation for the tunneling 
between the tip and the sample. Our results for conductance steps justify the analysis carried 
out by Hirjebedin et al. [Science 317, 1199 (2007)] of observed step-like conductances by inelastic 
electron tunneling from spin excitations in a single magnetic adatom using a simple spin matrix 
element. In addition, our result gives a simple expression for the magnitudes of conductance steps 
and their lateral spatial variation with respect to the tip position, which can be calculated directly 
from spin-polarized wave functions at the Fermi level of the sample. 
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A recent breakthrough in the field of atomic-scale magnetism of surfaces was made by 
the observation of inelastic electron tunneling (lET) of spin-excitations in single and chains 
of magnetic atoms adsorbed on a surface using a low teniperature scanning tunneling micro- 
scope (STM) in the presence of a high magnetic field jl|, 3, 1^. The spin-excitations show up 
as steps in the measured tunneling conductance with respect to the voltage across the tip 
and sample. An analysis of the observed conductance steps in terms of simple spin-model 
Hamiltonians has revealed detailed atomic-scale information about the magnetic properties 
of adsorbed magnetic atoms on surfaces and their magnetic interactions . This analysis 
showed also the presence of selection rules in lET. Furthermore, in the case of single mag- 
netic adatoms the relative magnitude of the conductance steps could be accounted for by a 

u 

simple spin matrix element taken from inelastic neutron scattering [3] . The detailed physical 
mechanisms behind lET from spin-excitations and what determines the magnitudes of the 
observed conductance steps are still unclear and need to be resolved. 

In this note we develop a simple single-electron theory for lET from localized spin ex- 
citations of a single magnetic adsorbate on a non-magnetic substrate based on the Tersoff- 
Hamann approximation for elastic tunneling 4, |5[ . The spin of the adsorbate is treated in 
the adiabatic approximation. The inelastic coupling between the tunneling electron and the 
spin is via the exchange interaction. This coupling is treated in a non-perturbative, impul- 
sive approximation that was first developed in inelastic scattering in nuclear physics |6| and 
later applied to rotational scattering in molecular physics [?! . We derive a simple expression 
for the conductance steps whose magnitudes are determined by the spin matrix elements 
used by Hirjibehedin j2| and the local density of states of spin-averaged and -polarized parts 
of the one-electron functions. 

In developing our theory for lET from spin excitations, we consider specifically a local- 
ized, collinear spin moment of a single adsorbed magnetic atom on a non-magnetic metal 
substrate. The electronic structure is modeled in an one-electron approximation as provided 
by density functional theory. We begin by characterizing the electronic states of the sample 
as a function of the spin direction with boundary conditions specified along a fixed quanti- 
zation direction. Next we introduce the Tersoff-Hamann approximation for the differential 
tunneling conductance between the sample with a fixed spin direction and a non-magnetic 
tip. Before generalizing the TH approximation to lET from the localized spin excitation we 
present the adiabatic approximation for the spin dynamics. We end by a short discussion of 
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the results. 

The locahzed spin moment of the magnetic adatom arises from an exchange interaction 
between the electrons. In discussing the form of the electronic states, we neglect the small 
interaction terms arising from any external magnetic field and spin-orbit interaction. The 
one-electron states of magnetic adatom with a spin moment along a direction h is then 
obtained by a spin rotation of the states when the spin is along the fixed spin quantization 
direction z. A spin moment distribution s{r)z oriented along z is obtained in an one-electron 
approximation by a spin-dependent one-electron potential, 

v{r) = voo{r)I + Avoir)! + Avs{r)a, (1) 
^ vo{r)I + Avs{r)a, (2) 

where az is a Pauli spin matrix. Here vooir)! is the spin- independent potential of the 
bare substrate and Avoir)! + Avs{r)(Jz is the locahzed adsorbate- induced potential. Far 
away from the adsorbate in the bulk region, the one-electron potential is spin- independent, 
v{r) ~ Voo(^)/ resulting in two-fold degenerate scattering states |'0(U,<t') with energy e^. These 
states are determined by the Lippmann-Schwinger equation using the adsorbate-induced 
potential, the one-electron Green function of the bare non-magnetic surface and incident 
spin-degenerate states in the bulk region of the substrate, 

(f,<7|0^,,,) = 5.,.'(^1V'f) , (3) 
The states are labeled by an index // and the spin index a' = ±1 and they are coUinear 

(^O-l^/x,^') = Sa,a'i^i,,a'{r)- (4) 

To ease the notation when discussing the form of the electronic states for an arbitrary 
direction of the spin moment, we pair the two degenerate spinor wave functions (r, cr 1-0^,0-') 
into 2x2 matrices as, 

^,^;aa'{r) ^ {f,(T\llj^c,a') (5) 

and these coUinearw ave functions can now be represented in the following compact form 

Mr^ = ^^{r-)+r;{f1 (6) 
= ipOu,{r)! + iJsu.{r)<7z ■ (7) 



3 



Here the incident wave function Eq. ([6]) is the compact form of the wave function 

in Eq.©, 

= V^;(r)/ . (8) 

and ip^^ir) is the scattered part of the wave function from the adsorbate- induced potential 
Avo(r)/ + AvsirJcTz- In Eq.(I71), i/^o^{r) and '^ps^^{f) are the spin-averaged and -dependent 
parts of the wave functions, defined as, 

^0M(r1 = ^(^,,.=i(r1 + ^,,.=-i(r-)) (9) 
i^s,{^ = ^(V'^,.=i(r-)-V^^,.=_i(f)) (10) 

In particular, the spin moment distribution s{r)z is determined by these wave functions as, 

s{r)=A J2 mU^^s.ir^} (11) 

Note that the relative phases of the wave functions il^^{r) in Eq. ffTOj) are fixed by the 
boundary condition in Eq. ([8]). 

Now we consider the spin-moment to be aligned along an arbitrary direction h and the 
corresponding one-electron potential ^(r; h) is simply obtained by a spin-rotation using the 
matrix U{n) that rotates a spin state \a = l)s = \z)s to In)^: \n)s = U{n)\z)s and 

v{f;h) = U{n)v{f)U^{n) (12) 
= vo{r)I + Avs{r)n.a, (13) 

where we in the last step have used U{n)azU\n) = n.a. In the notation introduced in 
Eq-®, the new stationary states |'0ai,o-(^)) in the potential v{f; n) that obeys the incoming 
boundary conditions set by the states in Eq. ^ are then obtained by the spin rotation, 

^^(f;n) = U{n)^^{r)U\n) (14) 
= ipQ^ir)! + ipsf,{r)n.a. (15) 

Now we are in position to obtain the tunneling conductance between a non-magnetic tip 
and the magnetic sample with a fixed spin direction fi using the TH approximation. This 
approximation is based on the Bardeen approximation for the tunneling amplitude and a 
s wave approximation for the electrons tunneling from the tip. The differential tunneling 
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conductance G{V) at small positive sample biases V and at zero-temperature is then given 
by fl, 

G{V) = ^ J2 \T^a'M^)Me^-eF-eV)6{e,-eF) (16) 

where ei? is the Fermi energy and the tunneling (transition) amplitude T^j_„i.y„{n) from a tip 
state |V'i/,(t) to a sample state \iljyy{n)) is given by 

Tt,a'-,ua{n) =ts{'ilj^,a'{n)\rQ,a){s\'4)J) . (17) 

Here is the spin- independent s partial- wave amplitude around the tip apex tq for a 

tip state \i^u,a) and tg is determined by the tip. Note that here should the sample scattering 
states satisfy outgoing boundary conditions instead of the incoming boundary conditions in 
Eq. ([3]). Inserting the tunneling amplitude in Eq. (IT7|) into Eq. (1181) . one recovers the TH 
result that the differential conductance is proportional to the spin-averaged local density of 
states (LDOS) at the tip apex as, 

G{V) = C{Tr[p{eF + eV- fo, h)] (18) 

where Tr is the trace over the spin components and 

p^>„{e]fQ,n) = ^(fo,cr'|'0^^",fi)(?/'^^//,n|fo,cr)(5(e - e^) (19) 

and all tip quantities appear in Ct defined as, 

a = ^|t.|'$^|(5|V'.)p5(e,-e^) . (20) 

u 

The spin indices are avoided explicitly using the matrix notation introduced in Eq. ^ and 
Eq. (fT9!) simplifies to 

p(e; fo, n) = ^ V^^(fo; n)^l{rQ] n)6{e - e^) . (21) 

The result for ^^(f; h) in Eq. fll4l) and the simple properties of the Pauli matrices give directly 
the spin-averaged LDOS as, 

Tr[p(6;f,n)] = 2j](|^o.(rlP + |V^.,(r-)p)5(e-e,). (22) 

Thus from Eq. fll8l) . one finds that differential tunneling conductance G{V) is independent 
of the spin-direction h as expected from spin-rotation invariance. Note that the tunneling 
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amplitudes in Eq. (fT7I) with respect to z is dependent on h . After presenting the model of 
the spin-dynamics, we will show that this dependence gives rise to spin-excitations. 

The localized spin of a single adsorbed magnetic atom on a non-magnetic metal substrate 
is treated in the adiabatic approximation. The motion of the spin moment distribution is 
assumed to be rigid and the electrons are assumed to follow instantaneously the motion 
of the spin direction h. Thus the ground state energy of the electrons for a fixed spin 
direction S = a/ S{S + l)n determines the effective Hamiltonian Tie// for the spin motion. 
A perturbative treatment of the weak spin-orbit interaction and the weak external magnetic 
field B gives in general rise to an Tie// that contains terms up to second order in S and B. 
The 25*+ 1 spin eigen states of Tie// with eigen-energies are here simply labeled by 
their increasingvalue with an integer n. For instance, in the case discussed by Hirjebedin 
and coworkers 2], this Hamiltonian has the form 

T^e// = gfisB ■ S + DSl + E{Sl - Si) (23) 

where g is the gyromagnetic ratio, /i^ is the Bohr magneton, D and E are magnetic 
anisotropy parameters arising from the spin-orbit interaction and the z axis corresponds 
now to the surface normal. Finally, note that the external magnetic field and the localized 
spin-orbit interaction that gives rise to the effective spin hamiltonian are very weak com- 
pared to ^Vs{r)a.h in Eq.(l2]) and their effects on the variation of the electronic states with 
the spin direction fi has been neglected. 

In many cases the variation of the electronic structure of the energy scale of the spin-wave 
excitations is negligible justifying an impulsive approximation for the lET. In analogy to 
the impulsive approximation for inelastic scattering j^, the transition probability from the 
initial state of an incident electron in the tip with the spin of the adsorbate being in its 
ground state |<l>o) to a final state of an outgoing electron in the sample with the spin of the 
adsorbate being in a state |$„) is obtained from matrix elements T^o-'n;i/(TO of the tunneling 
amplitude T^a';ua{'n-) in Eq. (fT7Il as, 

{'^n\T^a';ua{n)\%) (24) 

Taking into account the Pauli exclusion principle in the summation of tunneling probabilities 
over final electronic states, the resulting tunneling conductance will have a step AGn at a 
bias V corresponding to the spin-excitation energy AEn = E^ — Eq: eV = AEn. The 



magnitudes of these conductance steps are now obtained by replacing the matrix element in 
Eq. f[T6l) by the matrix element in Eq. fl24|) and are given by, 

= ^ \T^a'n;uao\'^S{eF - e^)5{eF - e^) (25) 

Note that the neglect of the small energy transfer and the use of the on-shell tunneling am- 
plitude is consistent with the impulse approximation. Using the TH result for the tunneling 
amplitude in Eq. (fT7|l and the matrix form of sample wave functions in Eq. (fT4|) . Eq. (!25|) 
simplifies to 

AGn = a$^Tr[($„|7/>^(fo;n)|$o)($o|V'i(ro;n)|$„)]5(ei. - e^) (26) 

AGn can now be evaluated using the explicit form for ip^{f\n) in Eq. fllSp . The matrix 
elements in Eq. fl2B]) are then obtained as, 

(<l>„|^,(fo; h)\^,) = ^o,{ro)Sno + ^.,(^0)^^^^^ (27) 

■y^b[b + 1) 

and from the elementary algebraic properties of the Pauli matrices, the trace of the matrix 
elements in Eq. (l26|) simplifies to 

Tr[($„|^,(fo;n)|<|.o)($o|^i(r-o;n)|$.)] = |^o.(ro) P^.o + \ts,{ro)\' ' ^fs + l)' ^^^^ 
which gives the following simple expression for the inelastic steps in the conductance. 

Here Po^s{^,ep) are LDOS, defined as, 

Po,.(r,e^) =2j2\As,{^MeF - e,). (30) 

The result for AGn in Eq. (l26l) is the main result of this paper and will now be discussed in 
some more detail. 

At biases larger than the maximum spin-excitation energy, the conductance reduces to 
the conductance Go for a fixed spin. In this bias range, the conductance is given by, 

2_^AGn = Gt \po[ro,tF) + ps{ro,tF) 2^ g(g+i) ) ^^^^ 



and using the identity, 

J2\{'^n\S\<^o)\' = S{S + l), (32) 

n 

which follows from the completeness relation, one obtains, 

Go = G{V = 0) = J2^Gn (33) 

n 

Thus the step conductancies AG„ normalized to Go are given by, 

^ = ^^''^ SiS + l) 
where the dimension-less strength parameter, 

70^0 = —j^ ^ p r, 35) 

only depends on the LDOS of the sample and is always less than unity. At low biases below 
the threshold for spin-excitations, the conductance (?< is given by AGq and 

Thus our results for AGn in Eqns flM|) and fl36l) justifies the analysis carried out by Hirjebedin 
et al. 0] of the observed step-like conductance by inelastic tunneling from the localized 
spin-excitations using | ($„|5'|$o) In addition, our result gives a simple expression for the 
magnitude of and its lateral spatial variation with respect to the tip position, which can 
be calculated directly from spin-polarized wave functions at the Fermi level of the sample. 

The result for lET from a localized spin excitation has some important similarities with 
the result from lET from a localized vibration. The step in the conductance AG when 
exiting a local vibrational mode of an adsorbate on a non-magnetic metal surface from its 
vibrational ground state to its first excited state was shown from perturbation theory and 



the TH approximation to be given by (8|, [s], 

AG = |(fo|5^^)|25(e^ - e,) (37) 

where is the first order change in the spin-degenerate one-electron states, |'0^(u)), of 



the sample with respect to a root-mean-square amplitude Su = y ^aM '^^ vibration with 
frequency fl, and mass M. This result for the conductance step can be written in a similar 
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form to the corresponding result for the spin excitation in Eq. (129|) but now using vibrational 
eigen-states as, 

AG = 2Ct 5^($i|(ro|^^(w))|$o)($o|(^M(«)ln))|$i)5(eF - e,) (38) 

Thus in this case the JET comes from a change of the amplitude of the tails of the wave 
functions in the vacuum region during the vibrational motion whereas in the case of a spin- 
excitation, the JET arises from the relative change in the amplitude between the two spin 
components along a fixed quantization direction of the vacuum tails of the spinor wave 
functions. 

Finally, we would like to make some comments on the underlying approximations and 
restrictions behind our result for the step-like conductance by inelastic tunneling from a 
localized spin excitation. The many electron effects on the lET are only introduced in 
the final state so that any influence of the opening of an inelastic channel on the elastic 
channel is neglected. In the case of lET from a vibrational excitation this influence can 
give rise to a decrease of the conductance instead of an increase |iq|. Furthermore, the 
spin-dynamics is treated in the adiabatic approximation so that life time effects resulting 
in a broadening of the conductance steps is neglected. The results can be straightforwardly 
generalized to non-zero temperature and negative bias. In analogous manner for inelastic 
electron tunneling from vibrations, the Fermi surface smearing of the tip and the sample 



give rise to a temperature dependent broadening of the steps with about bksT The 
results at negative bias is obtained by considering single-electron tunneling from the tip 
to the sample and gives that tunneling conductance G{—V) = G{V) is an even function 
with respect to the bias. The results can also be generalized to localized spin excitations 
of several adatoms and a non-magnetic tip but that and the application to specific systems 
are deferred to later work. 
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